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Diluted planar ferromagnets: nonlinear excitations on a non-simply connected 

manifold 
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We study the behavior of magnetic vortices on a two-dimensional support manifold being not 
simply connected. It is done by considering the continuum approach of the XY-model on a plane 
with two disks removed from it. We argue that an effective attractive interaction between the 
two disks may exist due to the presence of a vortex. The results can be applied to diluted planar 
ferromagnets with easy-plane anisotropy, where the disks can be seen as nonmagnetic impurities. 
Simulations are also used to test the predictions of the continuum limit. 
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Physics in two spatial dimensions has generated a lot 
of surprises. In condensed matter systems these surprises 
are specially interesting due to the possibility of techno- 
logical applications. The role of nonlinear excitations 
in the study of low-dimensional, artificially structured 
materials is a very exciting topic due to their observ- 
able effects on the physical properties of a realizable con- 
densed matter system. Nonconventional supports (such 
as curved surfaces) for condensed matter materials may 
induce a much richer physics because the interplay be- 
tween geometry and topology 1, 2, 3, 4]. Coulomb sys- 
tems living on a sphere |a,Jaj_magnetic systems living 
on a cylinder or a sphere 0, Qi Hi 13 etc have already 
been considered recently. Even flat surfaces (the plane 
R^) with pieces cut out from them open up new, inter- 
esting avenues of investigation. For example, if one or 
more sectors are excised from a single layer of graphite 
and the remainder is joined seamlessly, a cone results 
ISj ^J ■ By considering the symmetry of a graphite sheet 
and the Euler's theorem, it can be shown that only five 
types of cone can be made from a continuous sheet of 
graphite corresponding to the following values of cone 
angles 7 = 19.2°, 38.9°, 60.0°, 84.6°, 112.9°. These are 
all synthesized forming carbon nanocones Q- The coni- 
cal geometry is also interesting because of its connection 
with the problem of Einstein gravity in low dimensions 
[3. Another interesting possibility is the plane R^ with 
a disk of radius a cut out from it. Such a space is said 
to be not simply connected since it has the property that 
some closed curves drawn in it can not be continuously 
shrunk to a point. In two-dimensional magnetic materi- 
als described by the continuum limit, such a disk cut out 
from the plane has a simple physical interpretation. It 
can be seen as a nonmagnetic impurity present into the 
system [l^, [ill ll^ ll^ • Then, from now, it is convenient 
to use the language of magnetism. 



A large variety of two-dimensional magnetic materials 
are well described by the anisotropic Heisenberg models. 
Particularly, many of these materials are modelled as a 
continuum of classical spins with easy-plane symmetry. 
In these systems, vortices are important excitations, re- 
sponsible for interesting static and dynamic properties 
n^jTmEIH- In this letter, we explore the effect of the 
support manifold (the plane R^) not simply connected 
and focus on the vortex behavior. Specifically, we con- 
sider the plane R^ with two disks of radius a cut out 
from it. In the case of magnetism, the size a is the lat- 
tice constant. We note, however, that in experimental 
situations involving small scales (e.g., a small hole with 
size of the order of lattice spacing a) the applications 
of the continuum approximation (i.e., a long wavelength 
theory) should be viewed with caution. Then, we also 
compare some results of this theory with numerical sim- 
ulations on a discrete lattice. Our main motivation is to 
know how magnetic vortices could be pinned in the sys- 
tem. Recent works [l^, [HI [i2i US have shown that the 
topological excitations center is attracted and pinned by 
a nonmagnetic impurity. However, a real magnetic ma- 
terial contains many vacancies and then, the possibility 
of pinned vortices not centered on vacancies but among 
them should be investigated. Here, first we consider two 
spin vacancies and a single vortex to calculate the equilib- 
rium position for the vortex center. The configurations 
of energy minima are then obtained. It is shown that 
the effective potential experienced by a vortex due to the 
presence of the two lattice defects attracts the vortex cen- 
ter to either one of the vacancies center or the center of 
the line joining the two vacancies. Second, we argue that 
the presence of the magnetic vortex background may in- 
duce an effective interaction between the holes that is 
attractive. This interaction results in a geometrical frus- 
tration in the region between vacancies. 

The easy-plane ferromagnets are described by the 
Hamiltonian 
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(1) 



where J > is the exchange constant, < A < 1 is 



the easy-plane anisotropy and Si = {Sf,Sf,Sf} is the 
classical spin vector at site i. We focus on the case A = 
(XY-model). The spin field can be parametrized by two 
scalar fields 4> and m = cos9 (the azimuthal and polar 
angles) as follows S = {\/l — ni^cos(j>, \/l — ni^sincj), m}. 
The efi'ective interaction among a single vortex and 
the two spinless sites is discussed by using a simple 
continuum XY-Hamiltonian on a support manifold being 
not simply connected, written as 
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FIG. 1: A vortex (white circle) located at origin (0,0) near 
two holes (black circles) located at fi = {xi , 0) and r2 = 
(3^2, J/2)- The distance between the holes is p —\ f2 — ri \. 
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where the term in brackets is the Hamiltonian density 
H for the continuum version of pure Hamiltonian (1), 
the function Ui{r) defined as C/i(r) = 1 if |r — rij > a 
and Ui{r) = if |r — ri| < a, represents the spin va- 
cancy (the first disk) centered at position ri and the 
function U2{r) defined as C/2(r) = 1 if |r — r2| > a and 
U2{r) = if \f— r2\ < a, represents another spin vacancy 
(the second disk) centered at position r2. The distance 
between vacancies is p =| r2 — ^i \- We now make a 
very simple mathematical observation that follows from 
the associative property of the terms present in expres- 
sion (2) and that justify the equivalency between a hole 
on the plane and a vacancy in the Hamiltonian. Note 
that jH[Ui{r)U2ir)d'^r] = J[HUi{r)U2{r)]d^r and then, 
the problem of a magnetically coated plane not simply 
connected (left side of the identity) is equivalent to the 
problem of a magnetic plane simply connected but con- 
taining spin vacancies in the Hamiltonian (right side). It 
only depends on the point of view: the disks arc removed 
either from the plane or from the Hamiltonian density. 
With this in mind, we remember that simulations show 
that the vortex structure is not deformed appreciably due 
the presence of a hole [13 . It will also be assumed in our 
calculations. 

The XY-model supports only planar vortices [l9| . 
which the configuration is my = 0, (t>v = arctan[(y — 
yv)/{x — Xy)] for the vortex center localized at {xy,yy). 
The energy of this vortex configuration, in the absence 
of nonmagnetic impurities, is Ey — 7rJln(L/0.24a) |19| . 
where L is the system size. Then, considering the vortex 
center localized at origin (for simplicity), the effective in- 
teraction between the vortex and the pair of impurities 
is given by Veff{ri,r2) = E{fi,f2) - Ey, where E{fi,f2) 
is the vortex energy in the presence of the two impuri- 
ties placed at fi and r2 respectively. Using Hamiltonian 
(2) , the planar vortex energy in the presence of two holes 
(nonmagnetic impurities) is simply obtained leading to 
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where d is a suitable constant of the order of the lattice 
spacing a introduced in Eq.(3) in order to avoid spuri- 
ous divergences of the vortex energy (in the continuum 
approach) when the vortex center coincides with an im- 
purity center. Thus, to find this constant, we note that, 
in the case in which the centers of the three defects (two 
defects in the lattice and one in the spin field) are located 
at the same point (for example, at origin), then, the three 
" bodies" problem reduces to a simpler problem of a vor- 
tex on a single hole. In this case, 14// (0, 0) obtained from 
Eq.(3), should have the value -4.48 J (see Ref. [l3|) and 
hence, we get d = 1.1472a. Although the limit that the 
separation vanishes is reasonable in the continuum ap- 
proach, in a real discrete lattice two vacancies can not 
occupy the same point. Therefore, these approximations 
can give us important insights about the system with 
topological excitations and nonmagnetic impurities. 

The object now is to calculate the configurations of 
minimum energy involving these three defects. It means 
that we need to minimize the quantity 



Fd{ri,r2) = — 
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To do this, we consider the geometry shown in Fig.(l). 
In this figure, the vortex center is located at the origin 
and the coordinates of the nonmagnetic impurities are 
(xi, 0) and (a;2, 2/2)- Here we have chosen 2/1 = (i.e., the 
hole 1 located along the x-axis) for effect of simplicity. 
Then, Eq.(4) is rewritten as 
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Using the Lagrangian multipliers method, one has 

VFd(a;i,a;2,y2) = -KVip{xi,X2,y2), with ip{xi,x2,y2) = 
p2 — (a;2 — Xi)^ — y2 ~ 0, where (/s is a constraint that 
keeps invariable the distance p between the two vacancies 
and K is a parameter. Thus, we get the following system 



of equations to be solved for xi^ X2 and y2 



r 2 lA 2 2\ 



(a —d -Xi) I \ r\ 

^2 (^^■f+d2)ix-i+yi+d^)2 - K[X2 - a:i) = 

2 i2 2 

r a —d —x^ |_ 1 p. 

p2 - (x2 -xi)2 -2/1 = 

We found two solutions for the above equations that sat- 
isfy our requirement of extremuni for Fd : 
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{xi,X2,y2) == (0,0, ±p). 



(7) 



(8) 



The identification of these solutions as energy minima 
will follow from a consideration of the original equations. 
Substituting Eqs. (7) and (8) into Eq. (3), one obtains 
the two energies Ei^ and E2v as a function of p for these 
configurations, respectively 
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Hence, in general, any configuration of these three de- 
fects in which the vortex center is found at the central 
point of the line joining the two impurity centers, leads 
to the situation with energy Eiy. By the other side, con- 
figurations of this system in which the vortex center is 
found at the center of one of the two impurities, produce 
the energy E2v 

Now, it should be interesting to check the conditions 
of stability of these two possibihties. It can be done by 
comparing the energies Eiy and E2v as a function of the 
impurities separation p. Figure (2) shows the functions 
£,i{p), defined as 6(p) — Eiy — Ey, {i = 1,2). Note 
that ^it, > S,2v for any p > 1.087a and ^i„ < £,2v for 
< p < 1.087a. It indicates that configuration 2 (with 
the vortex center coinciding with an impurity center) is 
energetically favorable for an appreciable range of im- 
purities separation p. It is interesting to note that, in 
the case of the two nearest neighbor vacancies (separated 
by the distance a), the vortex center should occupy the 
center of the intersection of holes (or the center of the 
large plaquette containing the two neighbor spinless sites 
in a discrete lattice). Figure 2 also shows that the po- 
tential ^2 is practically constant for p > 2a, and hence 



£(J) 
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P(a) 



FIG. 2: Potentials ^i (dashed curve) and ^2 (solid curve) 
versus p. For p > 1.087a, ^2 < £,i- 



the force between the vortex-on-vacancy and the other 
vacancy decays rapidly with the distance of separation. 
Note also that, in this case, the energy necessary to re- 
move a vortex from a vacancy is 2.24J, which is smaller 
than 4.48 J, indicating that the presence of other impu- 
rities lower the biding energy of the vortex-on-vacancy 
state. Extrapolating the present results for quenched di- 
luted layered magnetic materials with very low impurity 
concentrations (in this case it is expect that the vortex 
density is almost the same as the vacancy density), it 
is conceivable that a magnetic vortex lattice (not neces- 
sarily periodic, since the vortices distribution would be 
alike the vacancies distribution) could be present into 
the system. Besides, this vortex lattice could not be so 
rigid, because, as we have seen, the impurity concentra- 
tion must decrease the vortex-on-vacancy binding energy, 
giving some mobility to vortices. This picture would have 
important consequences to the spin dynamics, which may 
be observed in experiments. For example, the neutron 
scattering function Sj-y{q,uj) may still exhibit a central 
peak (the cause of this peak is believed to be the vortex 
translational motion jl5^, 16, 17J), but its shape, width 
and height should modify considerably. 

Next, we investigate the nature of the interaction be- 
tween two static vacancies through the vortex back- 
ground. In Ref . |13l| . it was shown that the energy of 
a vortex at origin in the presence of only one impurity at 
distance | n \ away is given by Ey + (tt J/2) ln[l — a? /{r'j + 
iP')], where b — 1.03a. Hence, the energy cost of removing 
one spin (placed at ri) from an infinite plane containing 
a vortex at origin is ei(r'i) = — (7rJ/2) ln[l — a^/(r| + 6'^)]. 
Of course, due to the cylindrical symmetry, this energy 
does not depend on the direction of r^. From Eq. (3), it 
is easy to see that the energy cost of removing two spins 
(located at sites ri and r2) from the vortex background 
is 
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Thus, there is an effective interaction potential between 
two static holes through the vortex background given by 
A = e2(ri,r2) — ei(n.) — £1(^2)- This potential is, then, 
expressed as 



A = — In 
2 



l + a^fb(ri,r2) 
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where i^b(ri,r2) is also given by expression (4) but with 
d substituted by b. This effective interaction depends on 
the distance of separation as well as the orientation of 
vacancies in the vortex background. Since b < d, this 
energy is negative, which means that the effective inter- 
action between vacancies is attractive. The problem can 
be summarized as follows: consider an initial configu- 
ration of these three defects as shown, for example, in 
Fig. (1). Supposing that p > 1.087a, the spin distribu- 
tion (or the magnetic energy) will be then readjusted to 
minimize the system energy, making the vortex center to 
move to the nearest vacancy center. However, the spin 
distribution did not reach its minimum in energy yet. 
The vortex-on-vacancy is still attracted by the other va- 
cancy, but now, neither the vortex-on-vacancy state nor 
the other vacancy can dislocate on the plane. In fact, in a 
rigid plane (lattice), the holes are fixed to their original 
positions and can not move on the plane. As a conse- 
quence, the system cannot reach its minimum in energy 
anymore; there is a geometrical frustration, which stems 
from a constraint (an invariable distance p between the 
two disconnected holes). Of course, the static- vacancy 
problem is an important limiting model for understand- 
ing the interaction of mobile holes through the magnetic 
background. If the vacancies were mobile, an attractive 
interaction could result in holes pairing around the mag- 
netic plane. Probably, in an elastic support, a variation 
in the geometry of the plane, compensating for the con- 
straint, may lead to a lowering of the energy. It is ex- 
pected that if the rigidity condition were relaxed, an elas- 
tic Hamiltonian density introduced in the problem should 
stabilize the plane against arbitrary deformations. 

To test some results obtained by using the continuum 
limit, we have also performed spin dynamical simulations 
on a L = 20a square lattice to study the behavior of 
a single vortex with its center initially located at the 
center of the system (0,0), in the presence of two spin 
vacancies located at (a/2, — 3a/2) and (—a/2, 3a/2) re- 
spectively (see Fig. (3)). This configuration starts with 
the vortex localized exactly at the central point of the 
line joining the two vacancies. Note that the analyti- 
cal results indicate that it is an unstable situation, since 
p « 3.16a, which is bigger than 1.087a. Then, it is ex- 
pected that, in the simulations, the vortex center must 
move in direction to one of the spinless sites. To see 
this we have imposed diagonally antiperiodic boundary 
conditions 20] Sl+i.v = -Si,L-y+i, So,y = -SL,L-y+i 

and Sx,L+i — Sl^x+1,1, S^a — —Sl~x+i,l, foi' all 
1 < 2;, y < i, in order to keep the vortex structure stable. 
The discrete motion equation for each spin is given by |21| 
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FIG. 3: Initial configuration with a vortex located at ori- 
gin and two vacancies at (—a/2, 3a/2) and (ci/2, — 3a/2) re- 
spectively. The squares represent the spinless sites and the 
distance between them is p = 3.16a. The length scale is ex- 
pressed in units of a/2. Many other configurations were also 
studied confirming the analytical results. 
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FIG. 4: Configuration after 7 x 10 time steps. The vortex 
travelled from the central position between vacancies to one 
of the vacancies center. It is in agreement with analytical 
results. 



dSi/dt = Si X h, where h = —JJ2ii^j^x + S^f^y) and Cx 
and iy are unit vectors in the x and y directions, respec- 
tively. The motion equations were integrated numerically 
forward in time using a fourth-order Runge-Kutta scheme 
with a time step of 4 x 10~*J~^. We notice that after 
10^ time steps the position of the vortex center always 
reaches one of the impurity centers (Fig. (4)). On the 
other hand, if the holes are separated by one lattice spac- 
ing, the continuunr theory predicts that the vortex center 
should stay on the central point joining the holes. Our 
simulations also agree with this result as it can be seen 
in figures 5 and 6. Figure (5) shows an initial configura- 
tion with a vortex at origin and two neighbor vacancies at 
(3a/2, a/2) and (5a/2, a/2) respectively. After 10^ time 
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FIG. 5: Initial configuration with a vortex located at origin 
and two vacancies at (3a/2, a/2) and (5a/2, a/2) respectively. 
The distance between these two vacancies is la. 
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steps, the vortex center stays between the vacancies (Fig. 
(6)). However, there is the possibility that this last con- 
figuration may not be completely static, but the vortex 
center seems to oscillate around the central point joining 
the vacancies, affecting the spin dynamics. The analysis 
of this possibility is out of the scope of this work and will 
be investigated in a future paper. 

In summary, we have studied the problem of the mag- 
netic vortex behavior in a non-simply connected mani- 
fold. The results can be applied to magnetic materials, 
including antiferromagnctic systems. The flat plane con- 
sidered here was constructed by cutting two holes, which 
can be interpreted as nonmagnetic impurities. The cut 
edges (which are the boundary surfaces for the infinite 
two-dimensional plane) exert a profound influence on the 
vortex center. There are two ways of pinning a vortex, 
in which the stability depends on the distance p between 
the holes. Besides, we have shown that the vortex back- 
ground induces an effective interaction potential between 
two vacancies that is attractive. It generates geometri- 
cal frustration. If the holes and the resulting hole-vortex 
hybrids have metallic mobilities, then such an effective 
interaction could be a pairing mechanism for supercon- 
ductivity in doped layered antiferromagnets. The overall 
qualitative agreement between the analytical results us- 
ing the continuum approach on a plane with two holes 
and the simulations on a discrete lattice with two missing 
spins is striking. Our results may also have relevance to 
dilute two dimensional Josephson junction arrays. 
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FIG. 6: Final state of Fig. (5). Now, the vortex center stays 
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